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Abstract 

Distribution of fluorescence resonance energy transfer (FRET) efficiency between the 
two ends of a Lennard- Jones polymer chain both at equilibrium and during folding 
and unfolding has been calculated, for the first time, by Brownian dynamics simula- 
tions. The distribution of FRET efficiency becomes bimodal during folding of the 
extended state subsequent to a temperature quench, with the width of the distribution 
for the extended state broader than that for the folded state. The reverse process of 
unfolding subsequent to a upward temperature jump shows different characteristics. 
The distributions show significant viscosity dependence which can be tested against 
experiments. 
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1 Introduction 



Dynamics of polymer folding and unfolding in solution is a problem of much current 
interest. [|l], @, 0] Although some aspects of polymer folding appears to be under- 
stood, detailed experimental study of the folding scenario (especially the initial part) 
is still not available. Given the complexity of the problem, the computer simulations 
could consider only relatively simple models, such as a necklace of Lennard- Jones or 
square- well beads. In addition to its own intrinsic interest, the collapse of polymers 
from poor solvents has served as a theoretical model of protein folding in the early 
stages. I, I, |, 0, |, |, 0,0,0 

In a notable recent development, fluorescence resonance energy transfer (FRET) 
has been combined with single molecule spectroscopic (SMS) technique to provide 
a powerful novel approach to study the dynamics of folding. Deniz et al, |]T3 



re- 



ported studies of dynamics of protein folding by observing FRET in time domain 
from a single donor-acceptor (D-A) pair. In the measurement of Deniz et al, the 
suhpopulations |jl3|, |14|, [15], |l^ of the folded and denatured states of the protein 



chymotrypsin inhibitor 2 (CI2) were obtained as the concentration of the denatu- 
rant (guanidinium chloride) was varied. The most interesting result was that at the 
intermediate concentration of the denaturant, the distribution of Forster efficiency 
becomes himodal. It was concluded that this bimodal distribution is the signature of 
the "two state" nature of CI2 folding transition. We were curious to know whether 
such a bimodal distribution can be observed (a) for simple models like collapsing 
homopolymers and (b) whether it can originate from dynamics alone. To answer 
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these questions, we have carried out extensive Brownian dynamics (BD) simulations 
of FRET in single macromolecules. We find, surprisingly, that even such a simple 
system shows a bimodal distribution in the excitation transfer efficiency, for cho- 
sen sets of parameter values. The double peak is similar to the ones observed in 
experiments. [13] The bimodal distribution is observed not only during folding but 
also at equilibrium near the 6 temperature. The results obtained here should help 
in designing future experiments. 

The usually assumed mechanism for FRET is the Forster energy transfer (FET).|I8|, 



20 1 The rate of this transfer depends on the separation [R) between the energy donor 
and the energy acceptor. This rate, k{R)^ can be written as|T^, 

k{R) = k,{^^)' (1) 

where Rp is the Forster radius and kp is the radiative rate, which is typically in the 
range 10^ to 10^ sec~^ for the commonly used chromophores. The rate of energy 
transfer becomes equal to Kp when R = Rp. For a given D— A pair, Rp can be 
obtained by the usual method of overlap between the fiuorescence and the absorption 



spectra of the D-A pair.||T9[ For commonly used chromophores, Rp is fairly large, 
often as large as 50 A. This means that the rate is very large when the donor- 
acceptor pair is separated by a short distance. This may actually be a limitation of 
the Forster expression which is strictly valid when the separation between the D — A 
pair is much larger than their size. However, the above limitation shall have a minor 
effect in the present study which is qualitative in nature and aims at exploring the 
general aspects of the energy transfer efficiency distribution during folding. 
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In writing Eq.(l), the standard averaging over the orientations of the transition 
dipole moments has been carried out. In standard FRET experiments, the macro- 
molecule is doped with a D-A pair in suitable locations along the chain. Excita- 



tion transfer can be monitored by following the fluorescence either from the donor 
or the acceptor or from both. The time, Trxn, taken for the excitation transfer to 
occur depends strongly on the D-A separation R, as given by k{R) in Eq.(l). For a 
polymer (or a protein) in solution, both at equilibrium and during folding/unfolding, 
R is not only a fluctuating, stochastic function of time, but also varies in a definite 
way. In such cases, the distribution of the energy transfer efficiency contains non- 
trivial and useful information. Note that simulating FRET in real proteins is an 
exceedingly difficult problem - it is not trivial even for a simple homopolymer - a 
problem that has eluded theoretical description even today. 

We define the FRET efficiency {^f) by the following relation, 

$ ^(^) (2) 

k{R) + k 

rn 

where km is the rate which includes the radiative and non-radiative (other than 
Forster migration) decay rates of the donor-acceptor pair. We next define the FRET 
efficiency distribution P{^p) by the following expression 

P{^F) = ^rj:mF-^F{Tr.n)). (3) 
■'^ i=l 

The above equation is to be understood in the following fashion. After choosing 
a D-A pair at time t = 0, the pair is followed till the trajectory gets terminated 
due to the reaction. We define this time by Trxn- At this time, the existing end- 
to-end distance (i?) is used in Eq. (2) to calculate ^p. This was followed for J\f 



independent trajectories for pairs chosen from an equilibrium distribution. At the 
end we distribute the FRET efficiencies in to the bins of width 0.1. In this way, 
a continuous probabihty distribution {P{^f)) can be obtained from Eq. (3), by 
taking the A/"^ oo hmit - in our case, we get a histogram (figures 1-4). Similarly, 
we define the probability distribution of reaction times, 

1 ^ 

P{Trc,n) = T7 II ^(^ - ^rxn) (4) 
■'^ i=l 

In this paper, we present the calculations of the distribution of and Pijrxn) 

for Forster migration among polymer ends both at equilibrium and during folding 
and unfolding. In next section we describe the model and the simulation details. 
In section III we present the results. We close the paper with a few conclusions in 
section IV. 

2 Simulation details 
2.1 The model 

The model polymer chain is made of connected Lennard- Jones (LJ) beads. While 
this model homopolymer does not represent complex richness of a protein, it is 
known to show interesting folding kinetics. 0, |^ The total potential energy can be 
written as,p] 

TV i-1 N 

f^ = imwLj(r*i)+I]Mb(| -r,_i I) (5) 

i=2 j=l i=2 
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where N is the number of beads, rj is the position of bead i, Vij =| — rj \,ULj{r) 
is the LJ potential, 



ULj[r) = e 



(6) 



Ub is the bonding potential. 



-^ = ^E(|r.-r,|-6)^ (7) 

where a and e are the LJ collision-diameter and the well depth, respectively, fc^T 
is the thermal energy and k, represents the stiffness of the spring. Here, we use 
K, = 9, N = 80, and h = a. For convenience, we define e* = t/ksT. The unit of 
time, r, is / Dq. Thus, kp (= kph'^/Do) is also dimensionless. In this study we 
have chosen kp = 1 and 10. This choice of ^i? corresponds to the experimentally 
observed kp values, a bit biased towards higher values. For example, in a solvent 
with viscosity (r/) equal 1 cp, the radius of monomer molecules (i?) equal to 4^4°, 
kp = 1 corresponds to a kp of 0.76 ns while kp = 10 corresponds to that of 7.6 ns.||21|] 



R is scaled by b, the bead diameter, as usual. In viscous solvents, the viscosity can 
be much higher, and the kp can be even larger than 10. 

The time evaluation of the polymer chain is done according to the following equa- 
tion of motion, 

rj{t + At) = rj{t) + Fj{t)At + AX^{t) (8) 

where rj{t) is the position of j— th bead at time t. and the systematic force on j 
is denoted by Fj(t). The random Brownian displacement, AX*^(t), is taken from a 
Gaussian distribution with zero mean and 2 At variance. The time step. At, is varied 
from O.OOOlr to 0.0005r. The scheme of Noguchi and Yoshikawa0 has been used 

6 



to investigate the polymer folding. In this method, the equilibrium configuration 
obtained at e* = 0.1 is quenched by decreasing the temperature instantaneously to 
different values of e*, higher than 0.62 which is the 9 temperature in this model. 
The time dependent total energy, the root mean square end-to-end distance (i?^) 
and the radius of gyration (i?^) were all monitored to follow the progress of folding. 
The results presented here are the average over 10,000 of such trajectories with 
different initial configuration. More details on the simulation scheme can be found 



elsewhere. 1 22] To simulate FRET , we have probed many combinations of Rp and 
kp- We have selected Rp = 7, which is near the maximum in R^P{R) (we denote it 
by Rq), where P{R) is the end-to-end distribution. Another important parameter 
which affects FRET is kp. Large kp values result in the higher efficiency of FRET. 
In this study we have mostly dealt with kp=lQ. 

2.2 Time scales 



FRET in polymers involves several different time scales. Two time scales, Tpjoid 
and Tpurii are required to describe the average survival probability of FRET in 
the folded and unfolded states, respectively. ||22|] These two time scales are widely 
separated from each other, due to the sensitivity of survival time to the separation 
between the two ends. The third relevant time in this problem is TgfoM, the time 
required for the polymer to fold subsequent to the quench. For FRET to be useful 
in the study of folding, this Tqfoid should be intermediate and well-separated from 
the other two times. Two additional time scales kp'^ and b'^/Do come from Forster 
energy transfer and Brownian dynamics, respectively. While kp is fixed, B'^/Dq can 
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be varied by changing viscosity {rj). 



3 Results and Discussion 



Results for the distribution during folding process (subsequent to the quench from 
e* = 0.1 to e* = 0.8 ) are shown in figure 1. In this figure the simulated probabihty 
distribution of FRET efficiency during folding is plotted at Rp = 7 and 

kp — 10. One sees a clear bimodal distribution in the FRET efficiency. The first 
peak at low efficiencies arise from the extended state while the one at high efficiencies 
arise from the folded configurations. Note also that the distribution for the extended 
state is broad while that from the collapsed state is narrow. This is expected on 
physical grounds and has been observed in experiments. This bimodality is found 
to depend critically on the value of km which is a consequence of several competing 
time scales in the FRET. 

Results presented in figure 1 can be better understood from figure 2 where we have 
plotted the equilibrium FRET efficiency distribution in the extended (unfolded) and 
the collapsed (folded) states, at Rp — 7 and kp — 10. It is observed that the 
time taken for FRET is much larger in the unfolded state (Fig. 2a), compared to 
that in the folded state (Fig. 2b). This is reffected in the (ensemble averaged) 
survival probability {Sp{t)) of D-A pair (not shown), which is extremely fast in the 
collapsed state and very slow in the extended state. The reason for the observed 
(nearly) exponential distribution for the unfolded state (Fig. 2a) lies in the choice of 
Rp. Since maximum probability for end-to-end distance (obtained by maximizing 
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R'^P{R)) is at ~ 7.3, for the given temperature and interaction strength. Thus, 
at Rp — 7, there is significant population aheady at this distance. If we change Rp 
to smaU (hke Rp = 1) or large (like Rp = 10) values, the exponential distribution 
will be replaced by a Gaussian type distribution. This sensitivity of the distribution 
to Rp can be exploited in experiments. The distribution of reaction efficiencies in 
the folded state (Fig. 2b), however, is not exponential. The initial fast fall in the 
probability is followed by a somewhat slower decay. 

Both at high (e* 0.1 — 0.3) and low temperatures (e* ^ 0.8 — 0.9) the equilibrium 
P{^p) shows a single peak at lower and higher FRET efficiencies, respectively. How- 
ever, at intermediate temperatures (e* ^ 0.5 — 0.6), the equilibrium FRET efficiency 
distribution {P{^p)) again shows a bimodal distribution. This is shown in figure 3, 
where the equilibrium P($f) is plotted at e* = 0.5. The emergence of bimodality in 
P($f), even at equihbrium, is due to the closeness to the 9 temperature (e^ = 0.62). 
This indicates that a "two-state" model exists in this simple system of homopolymer 
chain, near the 9 temperature . 

The above observations seem to provide the following interpretation of the observed 
bimodality in figure 1. As the polymer collapses subsequent to the quench, the 
polymer passes through a succession of configurations. The initial configurations 
correspond to the extended state. The polymer passes through the intermediate 
state rather fast which shows the dearth of population at intermediate ^p, giving 
rise to the bimodality. 

The effect of viscosity (?]) can be studied by varying the dimensionless rate kp, 
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defined as kph^/Do- At constant kp, large values of kp represent solution of high 
viscosity and vice versa. In figure 4, P{^f) is plotted against the FRET efficiency 
at two very different values oikp. Open bars show the result iov kp = 1, while the 
filled bars represent that for the = 10. Figure 4 shows that viscosity can have a 
dramatic effect on FRET efficiency distribution. 

During the unfolding (when the temperature instantaneously raised to e* = 0.1 
from e* = 0.8), P($i;') shows a large peak at higher efficiencies which is accompanied 
by a tail (of smaller peaks) towards lower efficiencies. The probabihty distribution of 
reaction times ( subsequent to a temperature quench) P{Trxn) also shows a bimodal 
distribution. We have found that one can efficiently study the dynamics of the initial 
stages of folding by placing the pair near the chain end where the nucleation of the 
collapse starts in nearly all the cases that we studied. In this case distribution of the 
reaction efficiency is not as strongly bimodal as in the case previously mentioned. 

We have monitored the variation in the mean square radius and the energy during 
the polymer folding as a function of time. They show somewhat different behavior, 
(i?^) and (Rg) starts decaying only after an initial characteristic time delay. The 
total energy of polymer chain starts decaying immediately after the quench and 
continues to do so till it reaches the final stable minimum energy configuration. 
This is because initial decrease of energy does not require change in Rq - it occurs 
by establishing favorable contacts. 
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4 Conclusion 



To conclude, we have shown in this work that a bimodal distribution of excitation 
transfer efficiency and of reaction times emerge during folding and unfolding of model 
homopolymers in solution. The distribution looks surprisingly similar to the ones 
observed recently in the folding of real proteins by single molecule spectroscopy. The 
extent of bimodality is found to depend on the values of the Forster parameters {kp 
and RF),krn and the value of the diffusion coefficient (that is, the viscosity). Thus, 
a study like the one performed here can be useful in designing FRET experiments 
via single molecule spectroscopy. The present study suggests several exciting future 
problems. Both for folding and unfolding, one can initiate the FRET process after 
a suitable time lag (r) following the quench and can thus obtain a two dimensional 
distribution of P($f, t), like in NMR or ESR. FRET may also be able to differentiate 
between different collapsed states, like rod and toroid.0, |^ Simulations of this can 
be carried out by using a stiff polymer chain, p, |^ Simulation study by using more 
realistic models and also by incorporating the solvent molecules implicitly, may 
reveal more information. Further work in these directions is under progress. Finally, 
simulations of P{^f) during folding of real proteins remains a challenging task. 
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Figure captions: 

Figure 1. The distribution of FRET efficiency {P{^f)), subsequent to a temper- 
ature quench (from e* = 0.1 to e* = 0.8 ), obtained from BD simulations is shown 
at Rp — 7, kp —10 and N=80. This demonstrates the emergence of bimodahty in 
P{^f). Here km is fixed as 0.0001; the bimodahty is present at other values also, 
but sharpness depends on the magnitude of the radiative rate. 

Figure 2. The equilibrium FRET efficiency ($f) distribution is shown for (a) 
e* = 0.3 and (b) e* = 0.8. kp = 10 and Rp ^ 7 while km is fixed as 0.01. 

Figure 3. The equilibrium FRET efficiency distribution {P{(^p)) is shown for 
e* = 0.5, which is close to the 9 temperature (e^ = 0.62). kp, Rp and km are 
same as in figure 2. This figure demonstrates the emergence of bimodahty even at 
equilibrium. 

Figure 4. FRET efficiency distribution (P($ir)), subsequent to quench ( from 
e* = 0.1 to e* = 0.8) is shown for kp —1 (open bars) and kp — 10 (filled bars). 
Rp — 7 and km are same as in figure 2. 
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Fig. 1 Srinivas & Bagchi 




Fig. 2(a) Srinivas & Bagchi 




Fig. 2(b) Srinivas & Bagchi 




Fig. 3 Srinivas & Bagchi 




Fig. 4 Srinivas & Bagchi 
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